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Abstract 

Using the HERA - B detector, the bb production cross section has been measured in 
920 GcV proton collisions on carbon and titanium targets. The bh prodTiction was 
tagged via inclusive bottom quark decays into J/ip by exploiting the longitudinal 
separation of J/V' — > l^l~ decay vertices from the primary proton-nucleus inter- 
action. Both e'^e" and channels have been reconstructed and the combined 
analysis yields the cross section a{bb) = 32;ti2(stat) ^^(sys) nb/nucleon. 

Key words: hadroproduction, bottom mesons 
PACS: 13.85.Ni 13.85. Qk 13.20.He 24.85.-Fp 
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1 Introduction 



The theoretical description of heavy quark hadroproduction has been a subject 
of great attention in recent years [1^4]. For fixed target experiments [2,3], the 
predictions, based on perturbativc QCD, are compatible with experimental 
results from pion and proton beams, but both theoretical and experimental 
uncertainties are large. Only two measurements of the bb production cross 
section in proton-nucleus p + A interactions were previously reported [5,6]. 
Moreover, in collider experiments, the comparison between data and QCD 
predictions is not satisfactory [4]. 

The HERA - B experiment is designed to identify 5-meson decays in a dense 
hadronic environment, with a large geometrical coverage. Interactions are pro- 
duced on target wires in the halo of the 920 GeV HERA proton beam. The 
bb production cross section (cr^) on a nucleus of atomic number A is obtained 
from the inclusive reaction 

pA^bbX with bb ^ J/ijjY ^ (e+e-///+//-)y. (1) 

The 6-hadron decays into J/ ijj ( "6 ^ J/ il)" in the following) are distinguished 
from the large prompt J/%1) background by exploiting the b lifetime in a de- 
tached vertex analysis. We select b ^ J /ip decays in both the muon 
and electron channels and perform a combined bb production cross section 
measurement. 

In order to minimize the systematic errors related to detector and trigger 

efficiencies and to remove the dependence on the absolute luminosity deter- 
mination, the measurement is performed relative to the known prompt J/ip 
production cross section cxp [7,8]. Our measurement covers the J/ if] Feynman- 
X (xp) range —0.25 < xp < 0.15. Within our acceptance, the b to prompt cross 
section ratio can be expressed as: 

A^B _ Nb 1 

Aa^ Np ene^'Biibb^ J/i^X) ' 

where Aa^ and Aap are the b ^ J/ip and prompt J/ip cross sections limited 
to the mentioned xp range, Nb and Np are the observed number of detached 
b ^ J/ip and prompt J/ip decays. Er is the relative detection efficiency of 
b ^ J/ip with respect to prompt J/ip, including contributions from the trig- 
ger, the dilepton vertex and the J/ ip reconstruction, is the efficiency of the 
detached vertex selection. The branching ratio Br(66— > J/ipX) in hadropro- 
duction is assumed to be the same as that measured in Z decays, with the 
value 2- (1.16 ±0.10)% [9]. 
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Fig. 1. Side view of the HERA - B detector. 

The prompt J/ip production cross section per nucleon, a{pN J/ipX) ~ 
ap/A"", was previously measured by two fixed target experiments [7,8]. After 
correcting for the most recent measurement of the atomic number dependence 
(a = 0.955 ± 0.005 [10]) and rescahng [11] to the HERA - B c.m.s. energy, 
■s/s = 41.6 GeV, we obtain a reference prompt J/ip cross section of a{pN — >■ 
J/tjjX) = (357±8(stat) ±27(sys)) nb/nucleon. About 70% [7] of the J/i/j are 
produced in the kinematic range covered by our measurement. 

Since no nuclear suppression has been observed in D-meson production [12] 
and a similar behavior is expected in the b channel [13], we assume a = 1.0 
for the bb production cross section results presented here; i.e., = a{pN — > 
bb) ■ A. 



2 Detector, Trigger and Data sample 

A side view of the HERA - B [14,15] spectrometer is shown in Fig. 1. The 
spectrometer has a large geometrical coverage, from 15 mrad to 220 mrad in 
the bending plane and from 15 mrad to 160 mrad in the vertical plane. 

The target assembly [16] consists of two wire stations separated by 4 cm along 
the beam line, each containing 4 target wires of different materials. A servo 
system automatically steers the target wires during a run in order to maintain 
a constant interaction rate. The Vertex Detector System (VDS) [17] is realized 
by a system of 20 Roman pots containing seven planar stations (4 stereo views) 
of double-sided silicon micro-strip detectors (50 mm x 70 mm sensitive area, 
50 /xm pitch) that are operated in a vacuum vessel at 10 to 15 mm distance 
from the proton beam. An additional station is mounted immediately behind 
the 3 mm thick Al window of the vacuum vessel. 
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A dipole magnet with 2.13 Tm field-integral houses a first set of tracking 
stations, followed by a second set extending to 13 m downstream of the in- 
teraction region. To cope with the large particle flux gradient radial to the 
beam, the tracker is divided into a fine grained Inner Tracker (ITR) [18] and a 
large area Outer Tracker (OTR) [14,19]. The ITR uses micro-strip gas cham- 
bers (typical pitch of 300 /im) with gas electron multipliers. The OTR uses 
honeycomb drift cells with wire pitches of 5 mm near the beam and 10 mm 
away from the beam. 

The particle identification is performed by a ring imaging Cherenkov ho- 
doscope (RICH) [20], an electromagnetic calorimeter (ECAL) [21] and a muon 
detector (MUON) [22]. The RICH detector uses C4F10 as radiator. The fo- 
cal planes of the detector (above and below the beam line, respectively) are 
read out by multianode photomultipliers. The ECAL is based on "shashhk" 
sampling calorimeter technology, consisting of scintillator layers sandwiched 
between metal absorbers. In the radially innermost section, W is used in the 
absorber, and Pb everywhere else. The MUON system consists of 4 tracking 
stations located in the most downstream portion of the detector, at different 
depths in iron and iron-loaded absorbers. It is built from gas-pixel chambers 
in the radially innermost region and from proportional tube chambers, some 
with segmented cathodes (pads), everywhere else. 

The data sample presented in this analysis was acquired in a short physics 
run during the HERA - B commissioning period in summer 2000, at ~5 MHz 
interaction rate, with a maximum of two target wires operated simultaneously 
and separated by 4 cm along the beam direction. The two wires were made of 
carbon (1000 fira longitudinally and 100 /xm transversely) and titanium (500 
lim and 50 /xm, respectively). 

The data were collected by triggering on dimuon and dielectron signatures. 

The MUON pretrigger candidates were based on a double pad chamber co- 
incidence [23], while the ECAL pretrigger candidates were defined by ECAL 
clusters with a transverse energy Ex > 1.0 GeV [24]. The First Level Trig- 
ger (FLT) required two pretrigger candidates of the same type and forwarded 
these to the Second Level Trigger (SLT). The SLT is a software filter [25] run- 
ning on a farm of 240 PCs. Starting from the pad coincidences and high-E'^ 
ECAL clusters, a fast hit-counting algorithm and a simplified Kalman filter 
were applied to the OTR and VDS data to confirm the lepton pair candidates. 
An invariant mass cut of M > 2.0 GeV/c^ and an unlike- sign track require- 
ment were also applied in the electron channel. The data from accepted events 
were assembled and sent to the online reconstruction farm [26], consisting of 
100 dual-CPU PCs. The whole trigger chain allowed a reduction of the ini- 
tial interaction rate of 5 MHz to a final output rate of 20 Hz [27]. A total 
of ~450,000 dimuon and ~900,000 dielectron candidates were recorded under 
these conditions. 
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At the time of data taking, the ITR and the MUON pixel chambers were 
not included in the trigger. As a consequence, the forward hemisphere of the 
proton-nucleus c.m.s. is reduced in this measurement, compared to the full 
HERA - B acceptance. 



3 Monte Carlo Simulation 

A Monte Carlo (MC) simulation is used to determine the efficiency terms 
in Eq. 2 and to estimate the prompt J/ip background contribution to the 
b ^ J/t/j decay channel. 

The simulation of heavy quark (Q) production is achieved, first, by generating 
the basic process pA^" QQX including hadronization, using the Pythia 5.7 
event generator [28]; secondly, the remaining part of the process (X) is given 
as an input to the Fritiof 7.02 package [29] to simulate further interactions 
inside the nucleus. 

To describe the prompt J/ip kinematics accurately, the generated events 
are weighted according to the known prompt J/ip differential cross sections 
{da/dp^ and da/dxp) measured in proton-gold collisions [7]. These results 
were obtained in the positive Xp region, while our measurement covers the 
range —0.25 < xp < 0.15. MC studies based on the Color Octet Model [30] 
of charmonium production show a symmetric xp distribution of prompt J/ip 
decays. We therefore use the experimental parameterization [7] to extrapolate 
to the full xp space. The model dependence of the generated px spectrum is 
of less relevance since our pr acceptance is essentially flat. 

For the bb MC simulations, the events generated by Pythia are weighted ac- 
cording to a model with various contributions. First, the generated b quark 
kinematics {xp and px) are given by the computation of M. Mangano et al. [31] 
using the most recent next-to-next-to-lcading-logarithm (NNLL) MRST par- 
ton distribution functions [32] with a h quark mass of nib = 4.75 GeV/c^ and 

a QCD renormalization scale fj, — ^Jml -f- p^- Second, the intrinsic transverse 
momenta of the colliding quarks are smeared with a Gaussian distribution 
leading to {k^) = 0.5 GeV^/c^ [33]. Finally, the b fragmentation is described 
by a Peterson function [34] with a parameter e = 0.006 [5]. The subsequent 
6-hadron production and decay are controlled by the Pythia default parame- 
ters. The fo-hadron average hfetime is taken from Ref. [9]: ti, = 1.564± 0.014 ps. 

The sensitivity of the final result on the hb cross section within our acceptance 
(Act (66)) has been determined by varying the following hb MC model param- 
eters: the parton distribution functions (from MRST to CTEQ5 [35]), the b 
quark mass (in the range e [4.5,5.0] GeV/c^), the QCD renormalization 
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scale (from 0.5^7??.^+ to 2yml + pj^), the fragmentation function (from 
the Peterson form [5,36,37] with parameter e G [0.002,0.008], to the Kartvel- 
ishvih form [38] with parameter OjS = 13.7 ± 1.3 [37]), the intrinsic transverse 
momentum distribution (with (k^) in the range [0.125,2.0] GeV^/c^) and the 
fraction of 6-baryons produced in the b hadronization process in the range [0, 
12]%. The observed variations in the detection efficiencies have been included 
in the systematic error. 

The generated particles are propagated through the geometry and material 
description of the detector using the Geant 3.21 package [39]. A simulation 
of the detector response to particles is achieved by reproducing the digitization 

of electronic signals, with a realistic description of hit efficiencies and prob- 
lematic channels. The MC events are subjected to a full trigger simulation and 
reconstructed with the same algorithms as the data. 



4 J/iJ) Event Selection 

Since the observed number of prompt J/ijj decays is used as a normalization 
factor of the bb cross section measurement, we begin by selecting and counting 
the number of J/ip decays (iVp), before applying the detached vertex analysis. 
The lepton reconstruction in the OTR is seeded with the dilepton trigger track 
candidates; moreover a matching criteria is applied between the reconstructed 
track and the trigger track candidate in both the OTR and the VDS. The 
J/i/j ^ selection and counting procedure differs between the muon and 
electron channels, due to differences in the background levels, shapes and 
triggering conditions. 



Three criteria are used to select ^ decays and to purify the recon- 

structed sample from non-J/^ background: a dimuon vertex requirement and 
muon identification cuts in both the MUON and the RICH systems. The cuts 
are chosen to give the best signal significance {S/ on the number of seen 
J/ip (S) with the observed background {B). The resulting spectrum is shown 
in Fig. 2, with Np — 2880 ± 60 prompt J/ip ^ l^^pi''' decays. The like-sign 
spectrum shown in Fig. 2 is obtained from the same set of triggered events: 
the small discrepancy in number of reconstructed events in the background 
regions arises from the difference in trigger acceptance between the two cases 
and from physics contributions to the unlike-sign spectrum (Drell-Yan, open 
charm production). 
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Fig. 2. The /i"^//" invariant mass spectrum, after the J/ip selection cuts. The fit 
(soUd hne) assumes a Gaussian signal and an exponential background. The like-sign 
spectrum (dashed line) shows small discrepancies in the background regions (see 
text); it is not used in the analysis and serves only for illustration purposes. 

4.2 J/il) e+e- 

The selection of J/^ — > e'^e~ decays is more complex due to very large back- 
ground contributions, mainly from pions interacting in the ECAL and hadrons 
overlapping with energetic neutral showers. Due to such background, a clear 
J/i/j signal can be reconstructed only by means of strong electron identifica- 
tion requirements. Electron identification in the ECAL is based both on the 
ratio of the cluster energy E to the momentum p from tracking (E/p) and on 
the search for electron bremsstrahlung signals: 

- the E/p distribution is established for a purified J/ip ^ e^e^ sample by using 
a double-bremsstrahlung requirement as described in the following paragraph. 
The E/p spectrum is compatible with a Gaussian distribution of mean 1.00 
and width a ^9%; 

- bremsstrahlung photons emitted upstream of the magnet maintain the origi- 
nal electron direction; thus they can be used to correct the electron momentum 
at the vertex and they also provide a clean electron signature (bremsstrahlung 
tag). 

The e'^e~ invariant mass distribution is shown in Fig. 3(a), requiring only 
that E/p be within 1 a from unity for each track. Figs 3(b,c) show the im- 
provements in signal significance that are obtained when the bremsstrahlung 
selection is added to the E/p requirement. 

Table 1 lists the number of prompt J/ijj signal events found for different E/p 
and bremsstrahlung requirements. Using these J/ip sets one can measure the 
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>0 bremsstrahlung y a) 



>1 bremsstrahlung y b) 




2 bremsstrahlung y c) 



2 2.25 2.5 2.75 3 3.25 3.5 3.75 4 
e'^e" invariant mass (GeV/c^) 



4.25 4.5 



Fig. 3. The e+e^ invariant mass spectra, with an E/p cut at 1 a on both e+ 
and e~ tracks; (a) no bremsstrahlung requirement; (b) with at least 1 identified 
bremsstrahlung photon; (c) with two identified bremsstrahlung photons. The fits 
(solid lines) assume a Gaussian signal and a polynomial shape for the background. 
The background shapes (dashed lines) differ from the dimuon case (Fig. 2), due to 
Et and invariant mass cuts in the trigger. 



Bremsstrahlung 
requirement 


0.5(7 


E/p cut in 
la 


units of a 
2a 


3a 


None 


895 ± 75 


2553 ± 184 


5362 ± 292 


[5710 ± 380 ± 280] 


> 1 


519 ± 45 ± 11 


1420 ± 70 ± 48 


2851 ± 108 ± 48 


3304 ± 148 ± 188 


2 


106 ± 13 ± 2 


308 ± 24 ± 6 


587 ± 38 ± 14 


661 ± 40 ± 29 


^brems 


0.34 ± 0.03 ± 0.01 


0.35 ± 0.02 ± 0.01 


0.33 ±0.02 ±0.01 


0.33 ± 0.02 ± 0.02 



Table 1 



The number of prompt J/V' — e'^e~ with different bremsstrahlung requirements 
and different E/p cuts on both tracks. The bremsstrahlung tag probability for a 
single electron is reported for each E/p cut. The value at 3a in square brackets is 
obtained by extrapolation from samples with stronger electron identification cuts 
(see text). 
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bremsstrahlung tag probability for a single electron redundantly, resulting in 
an average Cbrcms = 0.34±0.02(stat)± 0.02(sys). This measured tag probability 
is in good agreement with expectations from MC simulations. The systematic 
uncertainties include the observed fluctuations when varying the fitting func- 
tions and range used to estimate the amount of J/%1) events in the invariant 
mass spectra. 

The good knowledge of the particle identification efficiencies allows us to infer 
the number of prompt J/%1) present in a sample where looser identification 
cuts have been applied and where no clear J/%1) signal is directly visible. Under 

such conditions, the normalization factor Np can be obtained while preserving 
reasonable statistics for the final detached vertex analysis, which relies only 
on the vertex separation cuts for the background rejection, as will be shown 
in Sect. 5. 

The total number of prompt J/ip in our sample with no bremsstrahlung tag 
requirement and with a Scr E/p cut is Np = 5710±380(stat)±280(sys) (entry 
in square brackets in Table 1). 



5 Detached Vertex Analysis 

The long decay length of 6-hadrons is used to separate the b ^ J/ip events from 
the prompt J/ip and to further reduce the non-J/ip background. The decay 
length (Az), defined as the distance along the beam axis between the J/t/j 
decay vertex and the closest wire (primary production point), is shown in Fig. 4 
for a purified sample of J/ip ^ e^e~ events (2 bremsstrahlung requirement). 
The width of the prompt J/iJj vertex distribution is more than 10 standard 
deviations smaller than the mean decay length of triggered 6-hadrons (r^O.S 
cm) . Given the achieved vertex resolution, a detached vertex cut proves to be 
efficient in the signal selection. Additionally, a cut on the minimum impact 
parameter of both leptons to the production vertex (/^) or to the wire (lyj) 
is applied in the detached vertex selection. The minimum impact parameter 
distribution of prompt J/ip ^ e'^e~ decay leptons (purified sample) to the wire 
is compared to the corresponding distribution from b ^ J/ip decay leptons 
(dashed line) in Fig. 5, illustrating the potential gain in 6 ^ J/ip signal- purity 
when applying a minimum impact parameter cut. 

The prompt J/ip which survive the detached vertex cuts cannot be distin- 
guished from b ^ J/ip events and their contribution to the detached b signal 
has to be determined from MC. A study has been performed to ensure that 
the simulation tool reliably reproduces the real data for the physical quantities 
defining the detached selection cuts, as illustrated by the good MC-data agree- 
ment in Fig. 4 and Fig. 5. In the simulation of the channel, a Gaussian 
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decay length 



(cm) 



Fig. 4. The decay length (Az) of purified J/ip ^ e+e^ events (2 bremsstrahlung 
sample, Ri 15% background) for real data and compared to MC prompt J/ijj events. 
The resolutions of both data and MC distributions are in agreement (715 =b 24 /xm 
and 721 it 10 ^m, respectively). The arrows mark the cuts applied in the detached 
vertex analysis. 




Data 

Prompt //\|/MC 
b^J/\]fXMC 



0.02 0.04 0.06 0.08 0.1 
Lepton impact parameter to the wire (cm) 

Fig. 5. The minimal impact to the wire (/^) of the two lepton tracks for purified 
J/ip ^ e~^e~ events (2 bremsstrahlung sample) data and in MC for J/ip ^ e~^e~ 
events. The expected shape from b ^ J/tp decays is shown in dashed line (arbitrary 
scale). The arrow marks the cut applied in the e'^e" detached event selection. 

distribution has been added to the standard MC track-slopes, increasing the 
slope errors by 20% in average, in order to match the parameters observed in 
data and resulting in a similar good MC-data agreement as in the e'^e~ case. 

The optimization of the detached vertex cuts is achieved by maximizing the 
ratio S/y/B, where S is the number of accepted signal events in the MC 
b ^ J/ip sample and B the number of background events in real data, observed 
in the whole upstream region and downstream in the side bins of the J/i/j 
invariant mass signal. The systematic errors in the final a{bb) result take into 
account the variations in the signal estimated under different sets of cuts found 
by the optimization procedure. 
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Upstream spectrum 



Downstream spectrum 




|l |l invariant mass (GeV/c ) 



Fig. 6. The upstream (a) and downstream (b) invariant /i^/U" mass spectrum after 
detached event selection. The downstream curve shows the result of the unbinned 
likelihood fit, in which the yields of background and signal contributions and the 
background slope, as defined in the text, were left as free parameters. 

The main background contributions expected in the detached sample are due 
to combinatorics (bad vertex or track reconstruction) and to double semilep- 
tonic cc or bb events. The combinatorial yield is estimated through the ob- 
served events in the region upstream of the primary interaction (unphysical 
region), while the charm background level is estimated by means of MC sim- 
ulations, assuming a cc production cross section of 40 /zb/nucleon for our 920 
GeV proton beam energy [40]. The bb background yield is estimated by study- 
ing the MC mass spectrum of bb events surviving trigger and selection cuts, 
with a contribution relative to the observed yield oib ^ J /ip events. 



5.1 



The detached vertex cuts found by the optimization procedure in the dimuon 
channel are: a minimum decay length of 7.5 times the uncertainty on the sec- 
ondary vertex position \^ a minimum track impact parameter to the assigned 
primary vertex of 160 /xm and a minimum track impact parameter to the 
assigned wire of 45 fim. 

Only 11 events survive these cuts downstream of the primary interaction region 
with invariant mass above 2.1 GeV/c^ (see Fig. 6). Only a single event is found 
upstream of the primary interaction. 

In order to use the full information of the detached invariant mass spectrum, 
an unbinned likelihood fit is performed, using the Gaussian parameters of the 
prompt J/ ip signal together with an exponential background contribution with 

^ The typical cutoff is at about 5 mm, with some events down to 3 mm being 
accepted. 
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free slope. The output of the fit shown in Fig. 6 yields l.OliJ h ^ J/il) events. 
The background slope obtained from the fit is compatible with the simulated 
charm and bottom quark background shape, although the statistics are low. 

Prom simulation, the expected background from prompt J/%1) decays is neg- 
ligible. The estimated background contributions of semilcptonic charm and 
bottom quark decays, together with the single event expected from combina- 
torial background (seen upstream) and the number of fitted signal events, are 
compatible with the 11 events observed downstream of the primary interaction 
region. 

To determine A(T{bb) in our xp range, the prompt J/ip and h ^ J/ip MC 
events are submitted to the same analysis chain used for real data. From 
simulation we obtain the efficiency terms entering in the cross section formula: 
Sr ■ e^^ — 0.41 ± 0.01. The corresponding bb cross section measured in the 
jji^ channel is A(t(66) = oj^jA = IG'^H nb/nucleon, obtained by using the 
weighted average of our target materials. All parameters contributing to the 
measurement (Eq. 2) are summarized in Table 2. 



5.2 b J/i/j e+e- 

The cut optimization procedure in the e"'"e~ channel results in the following 
criteria: a minimum decay length of 0.5 cm, a minimum track impact param- 
eter to the assigned wire of 200 //m or alternatively an isolation of the lepton 
candidate at the z of the wire from any other track by a minimum distance of 
250 /im. The types and values of the cuts are not exactly the same as in the 
muon analysis due to the very different background conditions. 

The detached selection yields 8 events upstream of the primary interaction re- 
gion (pure combinatorial background) and 19 downstream events (see Fig. 7). 
Among the downstream candidates, 10 events are found in the J/ijj mass win- 
dow (2.8 GeV/c^ < me+e- < 3.3 GeV/c^). 

Similarly to the muon analysis, an unbinned likelihood fit is performed on the 
invariant mass spectrum of the detached downstream e~^e~ candidates. The 
shape of the signal is taken from simulated b ^ J/ip decays, while the back- 
ground shape is a combination of the shapes obtained from simulated double 
semileptonic bottom quark decays and from pure combinatorial (upstream) 
events. The result of the likehhood fit is shown in Fig. 7(b), yielding 8.6l|;2 
b ^ J/ip events. When the background shape used in the fit is replaced by 
a pure combinatorial background shape or by a pure double semileptonic bb 
background, a ±7% variation is observed in the number of 6 — > J/ip events. 
This contribution is included in the systematic error of our measurement. 
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Upstream spectrum 



Downstream spectrum 




e'^e invariant mass (GeV/c^) 



Fig. 7. The upstream (a) and downstream (b) e+e" invariant mass spectra after 
detached event selection. The downstream curve shows the result of the unbinned 
likelihood fit, in which the yields of background and signal contributions, as defined 
in the text, were left as free parameters. 

The expected background from prompt J/'^ decays is of less than 0.2 events 
at the 90% C.L. in the whole downstream region. As in the fi^fi~ case, the 
estimated background yields from semileptonic charm and bottom quark de- 
cays, together with the expected combinatorial background level (8 events seen 
upstream) and the fitted signal, are compatible with the 19 events observed 
downstream of the primary interaction region. 

Different cut optimization techniques and assumptions have been tested to 
verify the stability of the signal. The optimizations are performed simultane- 
ously on the three detached vertex cuts {A.z, 1^ and the isolation cut) using 
the background from real data and the downstream b ^ J/ip ^ e^e~ events 
from MC. Independently on the optimization criteria, a J/ip signal with sig- 
nificance greater than 2 a is always observed in the downstream part of the 
spectrum, while a visible J/ip signal is never present in the upstream part. In 
Fig. 8, the selected detached events are displayed in a scatter plot of the in- 
variant mass versus the measured decay length (A^): a clustering is observed 
around the J/ip invariant mass for large A2; values in the region downstream 
of the primary interaction. 

To confirm the h assignment of the selected events, an unbinned maximum 
likelihood fit is performed on the decay length measurements. A mean decay 
length of 0.81 ± 0.03 cm on hb events is expected from MC. We measure 
1.0 ± 0.3 cm for the 10 downstream events in the J /ip region (2.8 GeV/c^ < 
rrie+e- < 3.3 GeV/c^), in good agreement with the bb interpretation, while 
the 8 upstream background events yield a mean decay length of 0.36 ± 0.13 
cm (measured using — A2;). To further verify that the selected events have 
features compatible with b decays, we performed a visual inspection of the 
candidates studying extra detached vertices (from the other b decay) and extra 
tracks attached to the J/V' vertex. Both categories of events are observed, and 
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Fig. 8. The scatter plot of e+e~ invariant masses versus the measured decay length 
(Az) for the selected detached events. The shaded region is removed by the Az cut. 
The horizontal line shows the mean J invariant mass value. A clear clustering of 
events around the J/V' mass with large A2; is observed in the downstream sample. 

their yields are compatible with MC expectations. Within the hmits of the 
available statistics, the J/ip xp and pt distributions are also compatible with 
the b ^ J/ip interpretation. 

Prom MC simulation, we obtain the efficiency terms entering in the cross sec- 
tion measurement (Eq. 2): er- e^^ — 0.44 ± 0.02. The corresponding bb cross 

section measured in the e~^e~ channel is Aa{bb) = a^/A = 38t\l nb/nucleon, 
obtained by using the weighted average of our target materials. All the pa- 
rameters used in Eq. 2 are summarized in Table 2. 



6 Combined Cross Section Measurement 

The two measurements, Aa{bb) — I6I12 nb/nucleon and Aa{bb) = 38li| 
nb/nucleon, obtained in the muon and electron channels, respectively, are 
compatible within statistical uncertainties. In order to extract the maximum 
information on the bb production cross section from our data, we combine 
the and e+e~ likelihoods in a four parameter likelihood maximization 

{Aa(bb), fi'^fi" background slope, and e~^e~ background yields) on the 

detached candidates. The fit provides our final result of the bb production cross 
section: 

Aa{bb) = 30li?(stat) nb/nucleon, (3) 

where the quoted uncertainty has been estimated directly from the fit (see 
Fig. 9). 
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1 9+?-? 8 6Hi 
-^•■^-1.5 "•'-'-3.2 


£r • £ n 


0.41 ±0.01 0.44 ±0.02 


Br{bb J/ipX) 


(2.32 ± 0.20 )% 


Aa{bb) 


16t\l iil>/iiucl. 38l|? iil-)/iiucl. 


Combined Act (66) 


30lif (stat) ± 6(sys) nb/nucleon 


Combined cr(66) 


32l}2(stat) ^^(sys) nb/nucleon 



Table 2 



The parameters entering into the o"(66) measurement (Eq. 2). 




bb cross section in our range (nb/nucleon) 



Fig. 9. The likelihood fits for the 66 production cross section in our xp range 
(Aa"(66)) using the fJ,'^fJ,~ and e'^e~ events separately (dotted and dashed line re- 
spectively) and in a combined analysis (solid line). 

The main sources of systematic uncertainty in the present measurement, which 
are not related to the final bb statistics, are due to the prompt J/ip cross section 
reference (11%), the branching ratio Br(66 — > J/ipX) (9%), the trigger and 
detector simulation (5%), the prompt J/V' MC production models (2.5%), 
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400 500 600 700 800 900 1000 

Proton energy (GeV) 

Fig. 10. The comparison of the HERA - B (2000) a{bb) value with other experiments 
and with the theoretical predictions of R. Bonciani et al. [2] updated with the NNLL 
parton distribution function in [32] (solid line: central value, dashed lines: upper and 
lower bounds) and N. Kidonakis et al. [3] (dot-dashed line: central value, dotted 
lines: upper and lower bounds). 

the bb MC production models (5%), the prompt J/tp —>■ e^e^ counting (5%) 
and the carbon-titanium difference in efficiencies (1.7%). Other contributions 
are below the 1% level. Uncertainties stemming from the background shapes 
used in the maximum likelihood fits on the invariant masses and from the 
cut values are dominated by the low statistics of observed detached events. 
For these sources we assign conservative uncertainties of ^24% ^ind 13% to the 
and e'^e~ channels, respectively. The overall systematic uncertainty for 
our measurement, averaged over the muon and electron channels, is of 

To compare our measurement with theoretical predictions, we extrapolate the 
Aa{bb) measurement to the full xp range, relying on the bb production and 
decay model described in Sect. 3 which foresee that 92% of J/ip from b decays 
are produced in our xp range. We obtain the total bb production cross section: 

a{bb) = 32+}^(stat) ^?(sys) nb/nucleon. (4) 

In Fig. 10, this result is compared with the latest QCD calculations [2,3] 
beyond next-to- leading order (NLO). The two predicted values at 920 GeV 
proton beam are, respectively, a{bb) = 25lill nb/nucleon and a{bb) = 
30 ± 13 nb/nucleon, in good agreement with our measurement. In the same 
figure, the E789 [5] and E771 [6] experimental results obtained with 800 GeV 
proton interactions on Au and Si, respectively, are plotted and are seen to be 
compatible (Fig. 10). 



^ Value based on Ref. [2], updated with the parton distribution function in Ref. [32]. 
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7 Conclusions 



Events coming from b ^ J/ip —>■ l^l^ decays liave been identified in a sample 
of ^1.35 million dilepton triggered events, acquired in a short physics run 
during the HERA - B commissioning period in summer 2000. The data analysis 
results in the identification of l.Ot^J 6 — J/ip ^ candidates and 8.6l'l3;2 

b ^ J/ip ^ e+e~ candidates. 

From these candidates, we compute the bb production cross section by nor- 
malizing to the known prompt J/ijj cross section. In the J/ip kinematic range 
-0.25 < xf < 0.15, we obtain Aa{bb) = IQtH nb/nucleon and Aa{bb) = 
38II5 nb/nucleon in the muon and electron channels, respectively. Within sta- 
tistical errors, the two results are compatible. The combined result of the bb 
production cross section measured by HERA - B at 920 GeV using pC and 
pTi interactions in our xp range is Aa{bb) = 30ln(stat)±6(sys) nb/nucleon. 
Extrapolating this measurement to the full xp range, we obtain the total bb 
production cross section: 

a{bb) = 32l}^(stat) ^^(sys) nb/nucleon. (5) 

This result is compatible with the existing measurements [5,6] and in agree- 
ment with the most recent QCD predictions [2,3] beyond NLO. 



Acknowledgments 

We express our gratitude to the DESY laboratory and to the DESY accel- 
erator group for their strong support since the conception of the HERA - B 
experiment. The HERA - B experiment would not have been possible with- 
out the enormous effort and commitment of our technical and administra- 
tive staff. It is not possible to list here the many individuals who have con- 
tributed to HERA - B. We are especially grateful to the following persons and 
their groups: G. Avoni, C. Baldanza, J. Bizzell, A. Cotta-Ramusino, F. Czem- 
pik, I. D'Antone, J. Davila, J. Dicke, A. Donat, U. Dretzler, A. Epifantsev, 
S. Fricke, W. Funk, A. Gutierrez, F. Hansen, M. Harris, S. Hennenberger, 
J. Hogenbirk, M. Jablonski, V. Kiva, M. Kolander, Y. Kolotaev, L. Laptin, 

H. Leich, H. Llidecke, Q. Li, K. Y. Liu, P. Liu, C. Lu, K. Ludwig, J. McGill, 
E. Michel, N. Murthy, E. Novikov, S. Omeltchuk, D. Padrazo, H. B. Pe- 
ters, P. Pietsch, M. Pohl, N. Ratnikova, A. Rausch, W. Reinsch, P. Rose, 

I. Rostovtseva, R. Rusnyak, W. Sands, P. Sole, S. Starostin^|], K.-H. Sulanke, 

^ deceased 



20 



V. Tchoudakov, M. Tkatch, K. Wagner, P. Wegner, V. Zerkin, E. Zimmer- 
Nixdorf. We thank the external HERA - B referees R. Forty, D. Froidevaux, 
R.-D. Heuer, K. Jakobs and J. Jaros for many stimulating discussions and 
suggestions. We are indebted to our administrative staff U. Djuanda and 
I. Kerkhoff for their continuous assistance. 

In the preparation of this paper, we have benefited from many useful discus- 
sions with M. Mangano, P. Nason, and R. Vogt on the theory of bb production. 



References 



[1] P. Nason, Proc. of the XX Int. Symp. on Lepton and Photon Interactions at 
High Energies, [hep-ph/01 11024 ; S. Prixione et al., Adv. Ser. Direct. High Energy 



Phys. 15 (1998), 609. 
[2] R. Bonciani et al., Nucl. Phys. B529 (1998) 424. 
[3] N. Kidonakis et al, Phys. Rev. D64 (2001) 114001. 

[4] S. Prixione et al., Proc of Int. Eur. Conf. on High Energy Physics, July 2001, 
Budapest, |hep-ph/01li368| . 



[5] D. M. Jansen et al., Phys. Rev. Lett. 74 (1995) 3118. 
[6] T. Alexopoulos et al., Phys. Rev. Lett. 82 (1999) 41. 
[7] M. H. Schub et al., Phys. Rev. D52 (1995) 1307. 
[8] T. Alexopoulos et al., Phys. Rev. D55 (1997) 3927. 
[9] D. E. Groom et al., Eur. Phys. J. C15 (2000) 1. 
[10] M. J. Leitch et al., Phys. Rev. Lett. 84 (2000) 3256. 
[11] T. Alexopoulos et al., Phys. Lett. B374 (1996) 271. 
[12] M. J. Leitch et al., Phys. Rev. Lett. 72 (1994) 2542. 



[13] R. Vogt, |hep-ph/0111271 



[14] E. Hartouni et al., HERA - B Design Report, DESY-PRC-95-01 (1995). 

[15] The HERA - B Collaboration, HERA - B Status Report, DESY-PRC-00-04 
(2000). 

[16] K. Ehret, Nucl. Instr. Methods A446 (2000) 190. 

[17] C. Bauer et al., Nucl. Instr. Methods A453 (2000) 103. 

[18] T. Zeuner, Nucl. Instr. Methods A446 (2000) 324; 
Y. Bagaturia et al., bep-ex/0204011, in press on Nucl. Instr. Methods A. 



21 



[19] M. Capeans, Nucl. Instr. Methods A446 (2000) 317. 

[20] J. Pyrlik, Nucl. Instr. Methods A446 (2000) 299; 
I. Arifio et al., Nucl. Instr. Methods A453 (2000) 289. 

[21] G. Avoni et al., Proc. of the IX Conf. on Calorimetry in Particle Physics, 
Annecy, France, Oct. 9-14, 2000, Calorimetry in High energy physics, (2001) 777; 
A. Zoccoh, Nucl. Instr. Methods A446 (2000) 246. 

[22] M. Buchlcr ct al., IEEE Trans, on Nucl. Sci. 46 (1999) 126; 
A. Arefiev et al., IEEE Trans, on Nucl. Sci. 48 (2001) 1059. 

[23] M. Bocker et al., IEEE Trans. Nucl. Sci. 48 (2001) 1270. 

[24] V. Alberico et al., II Nuovo Cimento, 110 A (1997) 1453; 
G. Avoni et al., Nucl. Instr. Methods A461 (2001) 332. 

[25] P. Kreuzer, Nucl. Instr. Methods A462 (2001) 212. 

[26] J. M. Hernandez, Proc. of the Int. Conf. on Computing in High Energy Physics 
CHEP 2001, Beijing, Sep. 2001. 

[27] J. M. Hernandez et al., Nucl. Science Symp. Conf. Record, 2000 IEEE, vol.2 
(2000) 12/116. 

[28] T. SjSstrand, Comp. Phys. Comm. 82 (1994) 74. 

[29] H. Pi, Comp. Phys. Comm. 71 (1992) 173. 

[30] P. Cho and A. Leibovich, Phys. Rev. D53 (1996) 150. 

[31] M. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B373 (1992) 295; P. Nason, 
S. Dawson and R. K. Ellis, Nucl. Phys. B327 (1988) 49. 

[32] A. D. Martin et al., Phys. Lett. B531 (2002) 216. 

[33] J. Chay, S. D. Ellis and W. J. Stirling, Phys. Rev. D45 (1992) 46. 

[34] C. Peterson et al, Phys. Rev. D27 (1983) 105. 

[35] H. L. Lai et al., Eur. Phys. J. C12 (2000) 375. 

[36] J. Chrin, Z. Phys. C36 (1987) 163; D. Decamp et al., Phys. Lett. B244 (1990) 
551; P. Nason and C. Oleari, Nucl. Phys. B565 (2000) 245. 

[37] A. Heister et al., Phys. Lett. B512 (2001) 30. 

[38] V. G. Kartvelishvih, A. K. Likhoded and V. A. Petrov, Phys. Lett. B78 (1978) 
615. 

[39] R. Brun et al., GEANT3, Internal Report CERN DD/EE/84-1, CERN, 1987. 
[40] S. Frixione et al, Nucl. Phys. B431 (1994) 453. 



22 



